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Abstract 

We find that the key features of the evolution and collapse of a trapped Bose 
condensate with negative scattering length are predetermined by the particle 
flux from the non-equilibrium above-condensate cloud to the condensate and 
by 3-body recombination of Bose-condensed atoms. The collapse, starting 
once the number of Bose-condensed atoms Nq reaches the critical value, ceases 
and turns to expansion when the density of the collapsing cloud becomes 
so high that the recombination losses dominate over attractive interparticle 
interaction. As a result, we obtain a sequence of collapses, each of them 
followed by dynamic oscillations of the condensate. In every collapse the 3- 
body recombination burns only a part of the Bose-condensed atoms, and Nq 
always remains finite. However, it can comparatively slowly decrease after 
the collapse, due to the transfer of the condensate paqrticles to the above- 
condensate cloud in the course of damping of the condensate oscillations. 

After the discovery of Bose-Einstein condensation (BEC) in trapped clouds of alkali 
atoms [0-0], one of the central questions in the field of Bose-condensed gases concerns 
the influence of interparticle interaction on the character of BEC. In this respect the Rice 
experiments with 7 Li attract a special interest, since a weakly interacting gas (n\a\ 3 <^1, 
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where n is the gas density, and a the scattering length) of 7 Li is characterized by attractive 
interaction between atoms (a < 0). As known [0, in spatially homogeneous Bose condensates 
with a < the negative sign and non-linear density dependence of the energy of interparticle 
interaction predetermine an absolute instability of the homogeneous density distribution, 
associated with the appearance of local collapses. A strong rise of density in the course of 
the collapse enhances intrinsic inelastic processes and leads to decay of the condensate. In 
a trapped gas the picture drastically changes. As has been revealed in [§,0, the discrete 
structure of the trap levels provides the existence of a metastable Bose condensate with 
a < 0, if the level spacing hu exceeds the interparticle interaction no\U\ (U = 4:7ih 2 a/m, m 
is the atom mass, n the condensate density, and u> the trap frequency). In terms of the 
number of Bose-condensed atoms N this condition can be written as 

N < iV 0c ~ l /\a\, (1) 

where Iq = (h/mco) 1 ^ 2 is the amplitude of zero point oscillations in the trap. 

The existence of a condensate in a trapped gas with a < has a clear physical nature. 
In an ideal trapped gas the BEC occurs in the ground state of the trapping potential. 
For attractive interparticle interaction the transfer of a condensate particle to exited states 
decreases the interaction energy by~no|t/|. But the energy of interaction with the trapping 
field increases by ~ fkv, and under condition ([!]) the change of the total energy is positive. 
In other words, there is a gap between the condensate and the lowest excited states. With 
increasing N to the critical value N 0c , the gap disappears and there will be an instability 
corresponding to the appearance of excitations with zero energy. As found in ||, for N < N Qc 
the non-linear Schrodinger equation for the condensate wavefunction \l/o has a stationary 
solution which becomes unstable at N >N 0c . 

However, the analysis performed in |7[] shows that the picture is more complicated. Ac- 
tually, for iVo < Nq c there are two global states with the same iVo and total energy E. In 
the first of them \l/o is almost a Gaussian and is localized in a spatial region of the size 
~ Iq. The other one is a non-stationary collapsing state localized in a much smaller spatial 
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region. The two states are separated by a large energy barrier and the transition ampli- 
tude is exponentially small, with the exponent depending on (N 0c — N ). From statistical 
considerations it is clear that in the course of accumulation of particles in the lowest trap 
level they turn out to be in the Gaussian state. There are peculiar fluctuations in this state, 
leading to the formation of "small dense clusters" of atoms. But the formation probability 
is again exponentially small (see 0). Thus, for Nq < Nq c the condensate will be formed in 
this metastable state which, however, does not decay on a time scale characteristic for the 
experiment. 

The problem of a metastable Bose condensate in a trapped gas with a < was also 
discussed in [pHTTH. The appearance of a Bose condensate with the number of particles 
N < N 0c was found in the Rice experiments with trapped 7 Li ||], where N 0c ~ 1000. 

The present paper is aimed at the analysis of the formation and evolution of a trapped 
condensate with a < in the presence of the Knudsen above-condensate cloud. Assuming the 
conditions which in the absence of collapse would provide the number of condensate particles 
N ^> N 0c , we show that the key features of the condensate evolution are predetermined 
by the particle flux from the above-condensate cloud to the condensate and by 3-body 
recombination of Bose-condensed atoms. Once the number of condensate particles reaches 
the critical value N 0c , the Bose-condensed cloud undergoes a collapse. However, we find 
that the compression reaches its maximum and turns to expansion when the density of 
the collapsing condensate becomes so high that the recombination losses dominate over 
attractive interparticle interaction. The recombination losses "burn" the condensate to 
Nq < N 0c , but the flux from the above-condensate cloud again increases N and a new 
collapse occurs et cet. As a result, we obtain a sequence of collapses, each of them followed 
by dynamic oscillations of the condensate. It is important that the recombination in the 
course of the collapse does not burn the condensate completely, and N always remains finite. 
As shown below, in a wide range of parameters the fraction of burned condensate particles is 
approximately one half. However, after every collapse N can further (comparatively slowly) 
decrease due to damping of the condensate oscillations, accompanied by the transfer of the 
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Bose-condensed atoms to the above-condensate cloud. 

We consider a Bose gas with a < and total number of particles N N 0c in an isotropic 
harmonic potential V(r) = muj 2 r 2 /2. The BEC transition temperature is determined by 
the relation T c = l.ObhcuN 1 ^ 3 [0 and greatly exceeds the interparticle interaction n\U\. 
Therefore, at temperatures T3>n| U\ (T <T C ) the equilibrium number of condensate particles 
can be found in the ideal gas approach: Nq = N[l — (T/T c ) 3 ]. The equilibrium BEC requires 
No smaller than Nq c , which immediately leads to the inequality AT = T C — T<T C . 

For iVo > Nqc the equilibrium BEC is impossible. Below we show that in this case there 
will be a strongly non-equilibrium evolving Bose-condensed phase. We discuss two limiting 
cases. First of them assumes that the conditions, which in the absence of interparticle 
interaction would lead to the equilibrium BEC with N ^> N 0c , are created abruptly. In 
this case, once the condensate is already present in the system, the flux of particles from 
the non-equilibrium above-condensate cloud to the condensate is induced by the condensate 
interaction with above-condensate atoms and is given by 

dN /dt = iN ; Y« 7o (l (2) 

where iV* is the total number of particles corresponding to N = N 0c . The parameter 7 is 
of order the frequency of elastic collisions and, hence, much smaller than u. The number 
of Bose-condensed atoms and, hence, the recombination losses per each collapse can not 
significantly exceed N 0c . Therefore, if initially N(t = 0) — iV* ^> N 0c , the major part of 
the condensate evolution proceeds with practically constant 7'. In the final stage, where 
N(t) — iV* is already comparable with A^ 0c , the decrease of 7' with N(t) becomes important 
and determines the approach of the system to the stationary regime. 

We will assume that the spherical symmetry of the Bose-condensed cloud, characteristic 
for iVo < Nq c , is retained when N reaches A^ 0c and the cloud collapses. A strong rise 
of density in the collapsing condensate enhances intrinsic inelastic processes. The most 
important will be the recombination in 3-body interatomic collisions. We will explicitly 
include this process (and the feeding of the condensate from the above-condensate cloud) 
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in the generalized non-linear Schrodinger equation for the condensate wavefunction. In the 
dimensionless form the equation reads 

i-^ = -A p ^o+p 2 ^o- |^o| 2 ^o-^o| 4 ^o+i7*o. (3) 

Here p = r/l , t = out/2 are the dimensionless coordinate and time variables, and — 
^o/n 1 / 2 , where the characteristic density h = (SnlUa])^ 1 ~ N 0c /8ttII. The recombination 
losses in the condensate and its feeding by the particle flux from the above-condensate cloud 
are described by the last two terms in the r.h.s. of Eq.(||]). The quantity £ = a r n 2 /uj, with 
a r being the rate constant of 3-body recombination, represents the ratio of the oscillation 
period in the trap to the characteristic recombination time 1 / a r n^ at no ~ n and the quantity 
7 = j'/u is the ratio of the oscillation period to the characteristic feeding time 1/7'. For any 
realistic parameters we have ( < 1. The parameter 7 is also small: As mentioned above, 
for the above-condensate cloud in the Knudsen regime one has 7' <C u. 

For N Q > Nq c Eq.(|3|) does not have stationary or quasistationary solutions even at 
£ = 7 = 0. Once the number of particles in the condensate exceeds the critical value N 0c , 
the Bose-condensed cloud starts to collapse. First it undergoes a purely dynamic compression 
determined by the non-linear interaction term — |\j/ | 2 ^o- The compression is accelerating 
with increasing \l/o> the compression time scale being r ~ l/|^ol- The total compression 
time is determined by a slow initial stage and turns out to be t ~ uj^ 1 (r ~ 1). From Eq.(||) 
one can see that the compression is constrained by the recombination losses and ceases when 
the condensate density reaches n ~ n* = U /ha r , i.e., 

i^oi 2 ~ i* *i 2 » r 1 > 1- (4) 

Three-body recombination accompanied by the particle losses predominantly occurs at 
maximum densities n~ra^. When the number of condensate particles becomes smaller than 
Nqc, the collapse turns to expansion and the 3-body recombination strongly decreases. The 
characteristic time interval, where the recombination losses are important, is ~ (a r n^) _1 
(t* ~ l^o*! -2 The total particle losses in the collapse are AN ~ N 0c a r nlt*. We see that 
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AN is independent of £ (and 7), or at least weakly depends on its value. This is confirmed 
by numerical calculations in a wide range of £. They show that the fraction of lost Bose- 
condensed atoms is approximately one half, although the internal structure of the collapse 
depends on the value of £. 

The recombination-induced turn of the collapse to expansion causes dynamic oscillations 
of the condensate: Due to the presence of the confining potential the expansion is followed by 
compression. These oscillations, with the period depending on u, resemble the condensate 
oscillations under variations of the trapping field (see fHI). We first perform the analysis, 
relying on Eq.(|3D and, hence, omitting the influence of the above-condensate cloud on the 
condensate oscillations. 

For revealing a qualitative picture we present the results of numerical calculation of 
Eq.(^|) with £ = 10~ 3 , 7 = 10™ 1 . Fig.l shows the time dependence of the number of Bose- 
condensed atoms, N (t). The time t = is chosen such that N (0) = 0.75N 0c and the 
Bose-condensed cloud is still stable with respect to collapse. The feeding of the condensate 
from the above-condensate cloud increases N Q and, once N becomes higher than N 0c , the 
collapse occurs. Three-body recombination in the course of the collapse burns approximately 
a half of the Bose-condensed atoms. Then, on a time scale ~7 _1 the particle flux from the 
above-condensate cloud increases the number of condensate atoms to iVo > iVo c and a new 
collapse occurs. It is accompanied by approximately the same particle losses as those in the 
previous collapse. The described oscillatory evolution of the condensate continues at larger 
times. The fine structure on the curve N (t), demonstrating moderate particle losses in the 
time intervals between the collapses, originates from the compression in the course of the 
dynamic oscillations of the condensate. 

The structure of the condensate oscillations is clearly seen in Fig. 2, where we present 
the spatial distribution of the condensate density, n (r, t)r 2 , at various times t. For t = t%, 
where the compression did not yet reach its maximum, the density at small r strongly 
increases compared to the initial distribution. But already after a short time (t 2 — ^1) "C uj^ 1 
the Bose-condensed cloud passes through the point of maximum compression, and both the 
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density and the number of condensate particles decrease due to recombination losses. Then 
the condensate starts to expand. A strong expansion of the condensate occurs at times of 
order u~ l (t = £3). The expansion is followed by compression, with a comparatively large 
increase of the density (t — £4). 

As already mentioned, the assumption of constant 7 relies on the inequality N(t)— ^> 
Nq c . When the latter violates because of the recombination losses, the parameter 7 decreases 
with N(t). In order to demonstrate the final stage of the evolution we present in Fig.3 the 
dependence N (t) calculated selfconsistently for the time-dependent 7, with 7' from Eq.(H) 
and N(t = 0)— iV* equal to 2.5iVoc and to 2Nq c . One can see that after two collapses the 
system approaches the equilibrium state, with N smaller than N 0c and depending on the 
value of N(t = 0) — N*. Again, N (t) always remains finite. 

One of the remarkable features of the collapse is the rise of dynamic energy of the 
condensate, induced by the recombination losses in the collapsing cloud. For No < N 0c 
the kinetic (K) and potential (P < 0) energy of the condensate are of the same order of 
magnitude, and the total energy E = (K + P) ~ N hu. In the course of the dynamic 
compression both K and \P\ strongly increase, whereas E is conserved. As a result, for a 
strong compression we have K, \P\ ^> E and, hence, K « \P\. Since K oc N , and \P\ oc Nq, 
the loss of 8Nq Nq particles changes the total energy by an amount 

6E=*^(2\P\-K)>0. (5) 

In the course of particle losses the relation between K and \P\ changes, which can reverse 
the sign of 5E. 

The process of 3-body recombination produces fast atoms and vibrationally excited 
molecules. Their kinetic energy is determined by the binding energy of the molecule and 
greatly exceeds the energy 5E/5N acquired by the condensate in the recombination event. 
Therefore, due to the recombination-induced increase of the condensate energy, the fast 
atoms and molecules simply carry away from the system slightly less energy than in the 
case of recombination in vacuum. Eq.(|]) and the above analysis implicitly assume that the 



mean free path of the fast atoms and molecules is much larger than the sample size and 
they escape from the trap without collisions with the gas atoms. The energy transferred to 
the system is concentrated in macroscopic oscillations of the condensate. In fact, this can 
be seen already in Fig. 2. 

It is worth noting that the recombination-induced increase of the condensate energy can 
lead to the appearance of short-wave excitations which overcome the trap barrier and carry 
away a significant part of the condensate dynamic energy. Together with a detailed analy- 
sis of damping of the condensate oscillations this problem is especially important for much 
smaller values of £ and requires a separate analysis. Here we only present qualitative argu- 
ments concerning the possibility of damping of the condensate oscillations. The damping is 
caused by the interaction of the oscillating condensate with the above-condensate cloud and 
for a large energy of the oscillations can be accompanied by the transfer of the condensate 
particles to this cloud. It will, certainly, occur on a time scale greatly exceeding the charac- 
teristic time of the collapse. On the other hand, the characteristic damping time is likely to 
be smaller than the time of feeding of the condensate from the above-condensate cloud (time 
interval between two collapses). Thus, effectively in each "collapse" the time dependence of 
the number of Bose-condensed atoms can consist of a sharp drop by approximately factor 
2, followed by a comparatively slow decrease of N . 

Let us now briefly discuss another limiting case, where the gas temperature is decreasing 
adiabatically slowly and for N (t) < N 0c the system is in quasiequilibrium. With decreasing 
T, the number of Bose-condensed atoms rises and, when it reaches N 0c , the collapse occurs. 
Similarly to the previous case, N drops. Since the total number of particles becomes smaller, 
the quasiequilibrium is reestablished at lower T c . As the temperature continues to decrease, 
N increases and the collapse occurs again et cet. This continues until N(t) > N 0c . It is 
important that for N(t) ^> iV 0c the instantaneous values of T and T c always remain very 
close to each other. 

Finally, we make a general remark. The collapse as a solution of the non-linear 
Schrodinger equation was a subject of extensive analytical and numerical studies. The 
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attention was focused on analyzing the character of the singularity and on finding univer- 
sal scaling solutions in the absence of dissipation or in the presence of weak dissipative 
processes (see [Q and references therein). Of particular interest was the search for the so- 
called strong collapse, which arrived at the concept of "burning point" (small spatial region 
absorbing particles). 

The picture of collapse, described in the present paper, stands beyond this analysis. To 
an essential extent this is related to the presence of the trapping potential which determines 
dynamical properties of the system and provides the existence of a peculiar quasistable con- 
densate with a limited number of particles. Another reason is that the collapse occurs in 
non-equilibrium conditions and there is particle exchange between the condensate and the 
above-condensed cloud. In other words, the condensate is an open system, which predeter- 
mines the appearance of a sequence of collapses. In this respect, BEC in ultra-cold trapped 
gases with a < opens possibilities for observing and studying novel pictures of collapse. 

We acknowledge discussions with R.G. Hulet. This work was supported by the Dutch 
Foundation FOM, by NWO (project NWO-047-003.036), by INTAS and by the Russian 
Foundation for Basic Studies. 
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FIGURES 

FIG. 1. The ratio N (t)/N 0c versus ut for f = 1(T 3 and 7 = HT 1 . The time t = is selected 
such that iV (0) = 0.75iVo c . 

FIG. 2. The condensate density profile for various times t. The dashed curve corresponds to 
t = (iV (0) = 0.75iV 0c ). 

FIG. 3. The ratio N (t)/N 0c versus ut for the time-dependent 7 (£ = 10~ 3 , 7(0) = 1CT 1 , 
A^o(0) = 0.75iVo c ). The solid curve corresponds to iV(0) — iV* = 2.5JV"o c , and the dashed curve to 
N(0)-N* = 2N 0c . 
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